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Resumen 


Canto-Ríos, J., Martín-Domínguez, A., Pérez-Castrejón, S., 
Alcocer-Yamanaka, V., € Tzatchkov-Velitchko, G. (mayo- 
junio, 2017). Modelación hidráulica de un reactor de 
electrocoagulación tubular de sección anular. Tecnología y 
Ciencias del Agua, 8(3), 111-126. 


La electrocoagulación es una opción en la remoción de 
contaminantes presentes en agua basada en el desgaste de 
electrodos para la formación in situ de un coagulante. En 
la eficiencia de este tipo de sistemas intervienen procesos 
fisicoquímicos, electroquímicos e hidráulicos. Diversos 
autores han estudiado esta tecnología desde los primeros dos 
enfoques, sin embargo, poco se conoce de la influencia de la 
hidráulica en la eficiencia del proceso. El presente trabajo se 
enfoca en el estudio de las pérdidas de carga totales en un 
reactor de electrocoagulación con flujo a presión y sección 
transversal anular simple y múltiple, donde esta última 
característica depende del arreglo de los electrodos. Se 
presenta el análisis de las pérdidas de carga por cortante 
y accesorios, las cuales se modelaron utilizando diferentes 
expresiones empíricas y se compararon con las mediciones 
realizadas a escala semipiloto. Como resultado del análisis se 
observó que la ecuación de Bahramir, Yovanovich y Culham 
(2006) asemeja las pérdidas reales en secciones anulares 
múltiples, mientras que la expresión propuesta por Davis 
describe mejor la sección anular simple. Estos resultados 
permitirán tener herramientas para diseñar reactores de 
electrocoagulación a escala real. 


Palabras clave: electrocoagulación, modelación, pérdidas de 
carga hidráulica, secciones anulares. 


Abstract 


Canto-Ríos, J., Martín-Domínguez, A., Pérez-Castrejón, S., Alcocer- 
Yamanaka, V., E Tzatchkov-Velitchko, G. (May-June, 2017). 
Hydraulic modeling of a tubular ring section electrocoagulation 
reactor. Water Technology and Sciences (in Spanish), 8(3), 111- 
126. 


Electrocoagulation is an option in the removal of contaminants 
present in water, based on the wear of electrodes for the in situ 
formation of a coagulant. In the efficiency of this type of systems are 
involved physicochemical, electrochemical and hydraulic processes. 
Several authors have studied this technology since the first two 
approaches, however, little is known of the influence of hydraulics 
on the efficiency of the process. The present work focuses on the 
study of the total charge losses in an electrocoagulation reactor with 
pressure flow and single and multiple annular cross section, where 
the latter characteristic depends on the electrode array. It is presented 
the analysis of the head losses by shear and accessories, which were 
modeled using different empirical expressions and were compared 
with the measurements made at semi-pilot scale. As a result of the 
analysis it was observed that the Bahramir, Yovanovich and Culham 
(2006) equation resembles the actual losses in multiple annular 
sections, whereas the expression proposed by Davis better describes 
the single annular section. These results will allow us to have tools 
to design real-scale electrocoagulation reactors. 


Keywords: Electrocoagulation, modeling hydraulic head loss, 
annular section. 
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Introducción 


La electrocoagulación es una alternativa viable 
para el tratamiento de agua, que se basa en la 
oxidación electroquímica de electrodos de sa- 
crificio para la producción de coagulantes que 
remueven contaminantes específicos. Tiene la 
ventaja de que no se agregan sales químicas que 
aumenten los sólidos disueltos (Holt €: Barton, 
1999), lo que reduce la producción de lodos de 
desecho. En la electrocoagulación se combinan 
tres procesos principales: fisicoquímicos, electro- 
químicos e hidráulicos (figura 1). La interacción 
de éstos ha dificultado el total entendimiento de 
los factores que impactan la efectividad de tal 
tecnología (Trompette $ Vergnes, 2008). 

El proceso hidráulico es fundamental en la 
eficiencia de los reactores de electrocoagulación 
debido a su influencia en las reacciones que se 
llevan a cabo en la capa límite de los electrodos 
y en la mezcla del coagulante con los contami- 
nantes presentes en el líquido. 

Los materiales para los electrodos de sacrifi- 
cio se seleccionan de acuerdo con la naturaleza y 
afinidad que tenga el contaminante que se desea 
remover con el coagulante producido (Piña- 


Soberanis et al., 2011); los de mayor uso son 
hierro (Fe) y aluminio (Al), y la configuración de 
los mismos puede ser plana o en forma de tubos 
concéntricos, dependiendo del tipo de reactor. 

Los reactores de electrocoagulación pueden 
clasificarse con base en sus configuraciones: 
cerrados o abiertos; y por el tipo de trayectoria 
que tiene el flujo: ascendentes, descendentes o 
de flujo mixto. Independientemente del tipo 
de rector, uno de los factores importantes para 
su diseño son las pérdidas de carga, las cuales 
son causadas por fricción o por piezas espe- 
ciales (Sotelo-Ávila, 2006). El presente trabajo 
revisa los modelos de pérdida de carga para 
seleccionar el más adecuado de acuerdo con la 
geometría estudiada. 


Antecedentes 


Hay pocos estudios del comportamiento hidráu- 
lico de reactores de electrocoagulación (REC) en 
la literatura. El cuadro 1 muestra algunas de las 
variables de interés que han sido estudiadas por 
diversos autores y que están relacionadas con 
el funcionamiento de este tipo de unidades de 
tratamiento de agua. 


Fisicoquímica 


Electroquímica 


Electro- 
coagulación 


Hidráulica 


Figura 1. Factores involucrados en el proceso de electrocoagulación. 
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Piña-Soberanis et al. (2011) revisaron 46 
estudios de electrocoagulación y llegaron a la 
conclusión de que no se describen con detalle 
los modelos experimentales empleados ni todas 
las variables que podrían ser útiles para tratar 
de interpretar los mecanismos en el sistema. 

Por otro lado, Khandegar y Saroha (2013) 
muestran algunas variables de operación para la 
remoción de contaminantes. Entre los resultados 
del análisis destaca la forma de los electrodos; 
concluyen que los electrodos con agujeros per- 
forados tienen mejor remoción que los planos 
y que estos últimos han sido poco estudiados. 
Mencionan que existe una separación óptima 
entre los electrodos, con lo que se tiene la mejor 
eficiencia; si la distancia es corta, los hidróxidos 
generados se degradan por colisión, si la distan- 
cia es mayor a la óptima, la eficiencia se reduce 
debido a la distancia que recorren los iones en 
el medio. 

Chen y Sheng (2004) mostraron que un área 
de contacto elevada entre los electrodos y el 
agua mejora la eficiencia de disolución de los 
mismos; sin embargo, las pérdidas de carga hi- 
dráulica aumentan debido a la fricción generada 
por las paredes internas del reactor. 

Martín-Domínguez, Rivera, Piña y Pérez 
(2008) evaluaron el impacto del gradiente de 
velocidad en el mezclado con respecto a la 
eficiencia de remoción de contaminantes en 
reactores de flujo pistón a presión y encontraron 
que la pérdida de carga aumenta debido a la 
acumulación de los hidróxidos en el interior del 
reactor, teniendo un efecto de taponamiento; 
dicho fenómeno se observa tanto en caudales 


bajos como elevados, disminuyendo la eficiencia 
del proceso. 

Entre las variables dependientes de estudio 
se encuentra la pasivación, la cual se refiere a la 
formación de una película relativamente inerte 
sobre la superficie de los electrodos, que afecta 
la reacción electroquímica. Asimismo, se tiene la 
densidad de corriente, que se refiere al flujo de 
electrodos que pasan a través de una superficie. 
Por último, se tiene la velocidad de mezclado, 
que se puede definir como la velocidad nece- 
saria para que dos componentes formen un 
producto homogéneo, que para el caso de la 
electrocoagulación se refiere a los flocs. 

Posteriormente, Piña-Soberanis et al. (2012) 
simularon hidráulicamente un reactor a pre- 
sión de sección rectangular utilizando dos 
paquetes comerciales basados en métodos 
numéricos avanzados para la simulación de 
problemas físicos, que resuelven las ecuaciones 
de Navier-Stokes en la simulación de fluidos 
(COMSOL-MultiphysicsO y PHOENICS-VRO) 
y compararon las pérdidas de carga hidráulica 
obtenidas con modelos matemáticos básicos 
(Dunlop-Nikuradse, Colebrook y Colebrook, 
Swamee-Jain); también analizaron las zonas 
muertas, y los cortocircuitos teóricos y reales. 

El objetivo de este estudio fue contribuir al 
conocimiento —en el ámbito de la modelación 
hidráulica— de reactores de electrocoagulación 
a presión, flujo pistón y sección transversal anu- 
lar simple o múltiple. 

Para lograr este objetivo se compararon di- 
ferentes ecuaciones utilizadas en el cálculo de 
pérdidas de carga, con la finalidad de encontrar 


Cuadro 1. Variables más estudiadas en REC. 


Variables independientes 


Variables dependientes 


Ph 
Flujo de agua a tratar 
Conductividad 
Orientación de los electrodos 
Disposición de electrodos 
Material de los electrodos 
Concentración del contaminante 


Corriente eléctrica 


Pasivación 
Eficiencia de remoción 
Densidad de corriente 
Velocidad de mezclado 

Voltaje 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 3, mayo-junio de 2017, pp. 111-126 


Canto-Ríos et al., Modelación hidráulica de un reactor 


de electrocoagulación tubular de sección anular 


aquella(s) que describiera(n) correctamente las 
pérdidas medidas en laboratorio para diferentes 
configuraciones de reactores. 


Marco teórico 


En los sistemas de conducción a presión, la 
energía total varía a lo largo de la trayectoria 
del flujo, lo anterior debido de manera principal 
a las pérdidas por cortante (/1) y las pérdidas 
locales provocadas por los accesorios utiliza- 
dos (h,) (Sotelo-Ávila, 2006). Las ecuaciones 
empleadas en este estudio consideran el factor 
de fricción de acuerdo con el material del ducto 
O la geometría de la sección transversal al flujo. 

Para el cálculo de la pérdida de carga en 
algunos modelos se requiere conocer el nú- 
mero de Reynolds, el cual se describe como la 
relación existente entre las fuerzas inerciales y 
las viscosas (ecuación (1)), también se necesita 
la rugosidad del material utilizado; esta última 
se puede consultar en tablas o por medio del 
diagrama de Moody (Sotelo-Ávila, 2006). Luego 
de comparar diferentes bibliografías, para este 
trabajo se consideró la rugosidad del aluminio 
como 0.015 mm y para el hierro 0.09 mm, como 
valores promedio: 


Re == (1) 


Donde v corresponde a velocidad del flujo; 
D, al diámetro de la tubería, y v a la viscosidad 
cinemática del agua. En el cuadro 2 se muestran 
los modelos utilizados en el presente trabajo 


para el cálculo de las pérdidas de carga por fric- 
ción. La ecuación (2) corresponde a la expresión 
desarrollada por Manning; la ecuación (3), a la 
de Hazen-Williams (Sotelo-Ávila, 2006), y la 
ecuación (4) a la de Scobey (Franqueti, 2005). 
Los modelos antes mencionados utilizan coe- 
ficientes de fricción fijos preestablecidos (CFF) 
y dependen sólo del material del cual se trate. 
La ecuación (5) corresponde a la expresión de 
Darcy-Weisbach (Sotelo-Ávila, 2006), cuyo coe- 
ficiente de fricción fconsidera características del 
material y la geometría de la sección transversal, 
que serán denominados coeficientes de fricción 
calculados (CFC). Finalmente, la ecuación (6) 
corresponde a la expresión de uso exclusivo 
para secciones anulares (Kartik, 2006). 

Cabe destacar que la expresión desarrollada 
por Darcy-Weisbach es la de mayor uso en el 
cálculo de pérdidas en ductos cerrados por las 
consideraciones hechas para su fácil aplicación; 
el factor de fricción utilizado en dicha ecuación 
se puede obtener con diferentes expresiones 
empíricas. El cuadro 3 muestra las expresiones 
utilizadas para calcular los factores de fricción 
para el modelo de Darcy-Weisbach en secciones 
circulares. 

Para el régimen laminar se utiliza la ecuación 
(7), que corresponde a la expresión desarrollada 
por Poiseuille en 1846 (Sotelo-Ávila, 2006); para 
el régimen de transición se usa la ecuación (8), 
que corresponde a la expresión desarrollada por 
Blasius en 1913 (Sotelo-Ávila, 2006); la ecuación 
(9) corresponde al factor de fricción propuesto 
por Nikuradse en 1920 (Sotelo-Ávila, 2006); 


Cuadro 2. Modelos de pérdida de carga. 


San (0?/D*>) 


(2) 


h,=10.674L[0') (Cy D*)] 


(8) 


h,=0.004098C,L (0 Ds) 


(4) 


(5) 


1 


Sa) (ym) (6) 
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Cuadro 3. Ecuaciones para calcular el factor de fricción en secciones circulares, modelo de Darcy-Weisbach. 


f=64/Re 


(7) 


f=0.3164 /Re"? 


(8) 


1/4/f =2108 (3.71D, /e) 


(9) 


1/4 --21og[[e/3.71 D,]+[2.51/Re |) (10) 


f=0.25/ (logíe /13.71D, +5.74./ Re"*]/2.302585)' (11) 


la ecuación (10) es la expresión de Colebrook- 
White de 1939 (Sotelo-Á vila, 2006), y finalmente, 
la ecuación (11) corresponde a la ecuación pre- 
sentada por Swamme y Jain en 1976 (Swamee 
éz Jain, 1976). 

Existen expresiones para el cálculo del 
factor de fricción de uso exclusivo para sec- 
ciones anulares (cuadro 4): Rothfus (1948), que 
corresponde a la ecuación (14); la ecuación (15) 
presenta la utilizada por Davis (Quarmby, 1966); 
la ecuación (16) muestra el desarrollo de Rehme 
(1973); los modelos de Muzychka y Yovanovich 
(1998), en la ecuación (17); Bahrami et al. (2006), 
con la ecuación (18); los parámetros importan- 
tes son el factor de forma (fr), ecuación (12); la 
relación de diámetros (k), en la ecuación (13): 


fría (12) 
d, 
a (13) 


En el cuadro 4 se encuentran las ecuaciones 
para el cálculo de los factores de fricción de las 
secciones anulares. 


Metodología 
Descripción 


La metodología para la obtención de datos de 
pérdidas de carga se dividió en tres partes: 1) 
pruebas en laboratorio con modelos físicos de 
sección anular simple; 2) cálculos en computa- 
dora con modelos conceptuales que utilizaron 
las ecuaciones descritas en el marco teórico; 3) 
aplicación del mejor o mejores modelos obteni- 
dos en los dos pasos anteriores, pero en uno de 
sección anular múltiple. 

Para las pruebas en laboratorio, primero se 
analizaron cuatro reactores de sección anular 
simple, R1, R2, R3 y R4, que variaron entre sí 
por el tamaño de los mismos y por el material de 


Cuadro 4. Ecuaciones para factor de fricción en secciones anulares. 


1/./f, = 4log(Re,/ ./f.)-04 (14) 
f=0.055Re""[(k -1)/k]” (15) 
V8/f=2.5InRe(£/8)+5.5-[F,, +2.51n2(1+K)] (16) 
fRe=8Vn|m(1+F?)/4E(e)] [1-57 (17) 

fRe =8 Vx (1-k)V1-k? /k?-1+(21n(£)+k?-1/1n(£)) (18) 
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los electrodos (aluminio o hierro). En cada uno 
de ellos se aumentó el gasto de forma paulatina, 
en un intervalo adecuado al tamaño de los reac- 
tores. Para el caso del R1, la variación del caudal 
para la medición de las pérdidas se realizó con 
un incremento de un litro por segundo; para R2, 
R3 y R4 cada 0.5 litros. Es importante recalcar 
que la tubería debía trabajar completamente 
llena, para asegurar la mayor área de contacto 
con los electrodos, por lo cual se tuvo un caudal 
mínimo de operación, como se muestra en el 
cuadro 6. Las pérdidas de carga entre la entrada 
y la salida de cada reactor se midieron con cada 
aumento del caudal, a fin de cuantificar las pér- 
didas de carga totales presentes en el modelo. La 
medición del caudal en cada prueba se llevó a 
cabo por el método volumen-tiempo a la salida 
del reactor en estudio, que consiste en dividir el 
volumen de agua recolectado en un recipiente 
entre el tiempo que transcurre en colectar dicho 
volumen. 

Para seleccionar el modelo o combinación de 
modelos de pérdida de carga que describieran 


el comportamiento de los reactores evaluados 
en laboratorio, se combinaron los modelos 
conceptuales, como se muestra en el cuadro 5, 
con el objetivo de utilizar cada uno de ellos en el 
intervalo de flujo para el que son válidos. En el 
cálculo de las pérdidas de carga y de los factores 
de fricción de acuerdo con el régimen de flujo y 
el modelo conceptual se programó un algoritmo 
en Visual Basic para Excel9, que combina los 
modelos descritos en el marco teórico. 

Los resultados de laboratorio se compara- 
ron con lo obtenido con los modelos teóricos, 
bajo las siguientes consideraciones: 1) el flujo 
a través del reactor fue de tipo pistón; 2) todos 
los elementos del fluido pasaban con la misma 
velocidad a través de todas las secciones (Váz- 
quez, 2005); 3) el flujo se consideró permanente 
y unidimensional, de manera que los perfiles 
de velocidad y la presión en cada sección de 
análisis fueran uniformes. 

Con la finalidad de seleccionar el modelo que 
mejor se ajuste a lo obtenido en el laboratorio, 
se compararon los resultados de cada expresión 


Cuadro 5. Modelos conceptuales utilizados para describir los regímenes de flujo. 


Modelo ID o 
Laminar Transición Turbulento 
1 MA Manning 
2 se Scobey 
HW Hazen-Williams 
Darcy-Weisbach (sección circular) 
P-D-N Poiseuille Dunlop Nikuradse 
P-C-N Poiseuille Colebrook-White Nikuradse 
C Colebrook-White 
P-C-S Poiseuille Colebrook-White Swamee éz Jain 
C-N Colebrook-White Nikuradse 
N Nikuradse 
ó P-D-B Poiseuille Dunlop Blasius 
Darcy-Weisbach (sección anular) 
P-M Poiseuille Muzychka $ Yovanovich 
P-BR Poiseuille Bahrami et al. 
R Rehme 
P-RO Poiseuille Rothfus 
P-DS Poiseuille Davis 
5 KJ Kartik 
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teórica con las pérdidas obtenidas en el modelo 
físico. Para lo anterior se obtuvo el error absolu- 
to (E), calculado como el valor absoluto de la 
diferencia entre la pérdida de carga teórica (),,) 
y la pérdida de carga medida en el laboratorio 
(h,), ecuación (19): 


abs 


E (19) 


abs r 


=(h,-5 


Dispositivos y materiales 


Los reactores utilizados en este estudio consis- 
tieron en módulos conformados por tubos de 
PVC unidos entre sí por codos de 90* y por un 
tramo recto de PVC (figura 2). En el interior de 
cada tubo se colocaron de manera concéntrica 
—para formar el polo negativo y el positivo— 
electrodos formados con tubos y varillas según 
el caso de estudio. Las configuraciones analiza- 
das fueron las siguientes: a) secciones anulares 


simples y b) secciones anulares múltiples (figura 
3). En ambos tipos de reactores se sellaron para 
cada reactor aquellas zonas donde no se quería 
que hubiera flujo de fluido. 

El cuadro 6 muestra las configuraciones y 
características de operación de los cuatro reac- 
tores (R1, R2, R3 y R4) de sección anular simple, 
que fueron utilizados para verificar la respuesta 
de los modelos estudiados en la primera etapa 
del estudio, así como el R5 de sección anular 
múltiple, con el cual se llevó a cabo la segunda 
parte del estudio. 

Para operar los modelos físicos se utilizó 
una bomba peristáltica con rango de trabajo de 
6 a 600 rpm, dando caudales de operación de 
16.80 a 1 680 ml/min utilizando tubería serie 
MasterFlex núm. 17 y una bomba sumergible 
para los reactores R1, R2, R3 y R4, así como una 
bomba centrífuga de % de pulgada para cau- 
dales de 30 hasta 250 1/min para el reactor R5. 
Las pérdidas de carga hidráulica se midieron de 


¡-— Módulo 


Tubo (PVC) 


Módulo 


Figura 2. Esquema de planta y corte longitudinal del reactor de electrocoagulación de sección simple. 
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a) b) 
Tubo 


Tubo sellado 


Tubo 
Carcasa 
(PVC) Carcasa 


(PVC) 


Varilla 


Sello 


Figura 3. Secciones transversales de las configuraciones estudiadas: a) anular simple y b) anular múltiple. 


Cuadro 6. Características geométricas de los reactores de sección anular simple. 


Sia Simple Múltiple 
R1 R2 R3 R4 R5 

Gasto (1/min) 10 - 35 15-95 1.0 - 6.4 0.9 - 7.8 30 - 207 
Número de Reynolds (adim) 2507-11384 | 1156-1 0391 770-8 382 709-10268 |3000- 34 000 
Material de los electrodos Fe Fe Al Al Aluminio 
Número de módulos conectados en serie 19 4 3 16 2 
Longitud del tramo recto entre codos (cm) 2.00 3.10 3.10 3.10 2.0 
Área efectiva de flujo (cm?) 7.08 1.30 0.88 0.88 13.03 
Perímetro mojado (cm) 29.74 6.85 5.70 5.70 61.69 
Radio hidráulico (cm) 0.24 0.19 0.15 0.15 0.21 

y] Diámetros (pulg) 

E Entrada y salida 11/4 3/8 3/8 3/8 2 

E Carcasa en PVC 21/2 3/4 3/4 3/4 21/4 

E Tubo 1 2 1/2 5/8 5/8 2 

a Tubo 2 11/4 

2 Tubo 3 3/4 

a Tubo 4 3/8 

E Longitudes (m) 

EE Carcasa en PVC 2.56 1.00 1.07 1.07 2.56 

E Tubo y varilla 3.00 1.07 1.12 1.12 3.00 

> 

E 

» la siguiente manera: en el R1, con manómetros Resultados 

3 de glicerina Instrutek, uno a la entrada de 2 

i kg/cm? y otro a la salida de 1 kg /cmí; en los La figura 4 muestra los números de Reynolds, 

3 reactores R2, R3 y R4 se utilizó un manómetro los cuales se calcularon considerando la tempe- 

S de mercurio. ratura del agua constante a 25 "C del R1 al R4, en 

E 
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Figura 4. Número de Reynolds para los diferentes casos de prueba. 


función de la velocidad de flujo correspondiente 
a los gastos de estudio, y con respecto al tipo de 
material de los electrodos de cada reactor. En 
la gráfica se observa que gran parte de la expe- 
rimentación se realizó entre los regímenes de 
transición (2 300 < Re < 4 000) y turbulento (Re 
> 4000), y que la pendiente de los dos reactores 
con electrodos de hierro es mayor que la de los 
de aluminio, debido a que los diámetros hidráu- 
licos de los reactores de hierro son mayores que 
los reactores de aluminio. 

La figura 5 muestra los números de Reynolds 
con respecto a la velocidad teórica de las cuatro 
zonas de flujo en el R5, y los reactores R3 y R4 
que tienen electrodos de aluminio; se observa 
que la zona tres tiene un comportamiento simi- 
lar a los primeros reactores estudiados debido 
a que el tamaño de la zona es similar a la del 
R3 y R4. 

En la figura 6 se compara el comportamiento 
del número de Reynolds con respecto a la ve- 
locidad del flujo de los reactores R1 y R2 (con 
electrodos de hierro) respecto al R5. En ella se 
muestra que la segunda zona y la zona central 


de flujo del R5 tienen similitud hasta el flujo 
máximo de operación de los primeros, pues su 
área efectiva de flujo es similar en dichas zonas. 

Con relación a la aplicación de los modelos 
conceptuales, la figura 7 presenta el promedio 
de error absoluto calculado para los cuatro reac- 
tores de sección anular simple, en función de los 
modelos seleccionados. Es posible observar que 
la combinación del modelo P-RO presenta el 
máximo error de ajuste con respecto a los datos 
de laboratorio. Las combinaciones P-D-N, P-C- 
N, E-N y sólo N presentan comportamientos 
similares y son los de menor error. 

Una vez analizados los errores absolutos 
para cada combinación, se revisan de manera 
independiente los errores absolutos presentes 
por reactor (cuadro 7). Se puede observar que 
el reactor con menor error es el R4; los modelos 
con menor error son P-M, P-BR y P-DS, con 
errores de 0.1 para los primeros dos y de 0.03 
para el último. 

Para un análisis más detallado de los errores 
que se obtienen con los modelos, al comparar 
el promedio de los errores absolutos (E_, ) de 
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Figura 5. Número de Reynolds por zona de R5 comparado con el de aluminio. 
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Figura 6. Número de Reynolds por zona del R5 comparado con los de hierro. 


los cuatro reactores de sección anular simple, en (figura 8); por lo tanto, independientemente del 
función del régimen de flujo (figura 8), se obtu- modelo de pérdida utilizado, los errores fueron 
vo lo siguiente: en laminar 0.017 m y en transi- menores en flujo laminar y de transición, por lo 
ción de 0.04 m, ambos con una varianza del 1%; que el modelo para el flujo turbulento es el que 
en turbulento, 0.146 m con el 28% de varianza define la combinación a elegir. 
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Figura 7. Promedio de error absoluto con los diferentes modelos para los reactores de flujo anular simple. 


Cuadro 7. Promedio de errores absolutos de pérdida de carga, por reactor y combinación de modelos. 


Rear R1 R2 R3 R4 
Modelo (m) 

P-D-N 0.58 1.05 1.35 0.19 
P-C-N 0.62 1.06 1.36 0.19 
C 0.43 1.12 1.47 0.34 
P-C-S 0.47 1.13 1.49 0.34 

C-N 0.62 1.06 1.36 0.19 
N 0.54 1.05 1.35 0.18 

P-D-B 0.78 0.99 1.23 0.31 

P-M 3.14 0.65 0.6 0.1 

P-BR 3.14 0.65 0.6 0.1 

R 3.24 1.42 205 0.82 

P-RO 12.21 2.66 4.31 2.29 

P-DS 2.66 (017/22 0.73 0.03 


La figura 9 muestra los promedios de erro- en todos los regímenes de flujo, los modelos 
res absolutos por tipo de modelo para cada para sección circular son los que tienen mejor 
uno de los reactores en estudio de acuerdo con ajuste. Como se observó en la figura 8, el ré- 
el modelo de cálculo (circular o anular), y se gimen turbulento presenta mayor error, que 
puede observar que al considerar las pérdidas afecta el promedio general de los modelos; 
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Figura 8. Promedio de los errores absolutos entre las pérdidas por fricción teóricas y reales, en función del régimen de flujo, 
para los reactores de flujo anular simple. 
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Figura 9. Promedio de errores absolutos en función del modelo utilizado y el reactor de estudio. 
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así, el análisis debe incluir los regímenes 
de flujo. 

Las pérdidas de energía calculadas teórica- 
mente resultaron ser superiores o inferiores a las 
obtenidas en laboratorio, principalmente para 
flujo turbulento; la selección de la rugosidad 
en las diferentes expresiones para el cálculo del 
factor de fricción puede derivar en la sobrees- 
timación o subestimación de la pérdida. Para 
establecer esto se compararon los modelos para 
el cálculo de la pérdida de carga que utilizan 
un factor de fricción preestablecido (CFF) y los 
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modelos que lo calculan (CEC). En la figura 10 
se muestra el resultado de esta comparación, 
considerando el tipo de factor de fricción y el 
material del electrodo. 

De lo anterior es posible observar que para 
cualquiera de los modelos y factor de fricción 
utilizado, los menores errores se obtienen con 
los electrodos de aluminio (menor área efectiva 
de flujo) y la mejor opción para este material 
son los modelos de sección anular con factor de 
fricción calculado. En el caso de los electrodos 
de hierro (mayor área efectiva de flujo), son 
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Figura 10. Promedio de errores absolutos en función del factor de fricción, modelo y material, para los reactores 


de flujo anular simple. 


los modelos para sección circular con factor 
de fricción calculado los que mejor resultados 
dan. Con base en lo anterior, lo modelos para 
sección circular (P, D, C y N) calculan de manera 
adecuada la pérdida de carga. 

Una vez seleccionados los mejores mode- 
los, se probaron en un reactor de aluminio de 
secciones anulares múltiples, el cual tiene las 
características mostradas en el cuadro 6, cuya 


geometría corresponde a la figura (b) de la 
figura 3. 

La figura 11 muestra los errores absolutos 
presentes para cada uno de las combinaciones 
de modelos propuestos. Se puede observar que, 
para este caso, los modelos para sección anular 
presentan mejor aproximación en relación con 
los circulares. Las combinaciones P-M, P-BR 
y HW son las que presentan el menor error, 


Promedio de errores absolutos (m) 


Modelos de pérdida de carga 


Figura 11. Promedio de errores absolutos en el reactor R5. 
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seguidas por P-DS y KJ, ambos grupos con 
errores menores a 0.50 m. 

La primera zona de flujo en el reactor es la 
comprendida entre la carcasa de PVC y el tubo 
1; la segunda, la formada entre el tubo 1 y el tu- 
bo 2, y así consecutivamente hasta el tubo 4. Es 
importante señalar que por los gastos utilizados 
en este reactor, los flujos son predominantemen- 
te turbulentos en todas las zonas del mismo. 

Debido a que los modelos anulares tuvieron 
un resultado con mejor ajuste a lo medido en 
laboratorio, se procede a aplicar las distintas 
combinaciones de modelos en cada una de las 
zonas de flujo, de manera que se seleccione 
aquella que resulte con menor error absoluto. 

De acuerdo con los resultados mostrados en 
la figura 12, la expresión que mejor se ajusta a 
lo medido en laboratorio es la de Bahrami et al. 
(BR), con un promedio de errores absolutos de 
0.790; seguida de Muzychka y Yovanovich (M), 
con 0.965; Hazen-Williams (HW) con 2.363, y 
Kartik (KJ) con 4.020 m. 

El aumento de las secciones anulares al inte- 
rior del tubo incrementa los radios hidráulicos y 


reduce las secciones de flujo, aumenta la fricción 
asociada con la pared de flujo y la rugosidad 
se incrementa notoriamente; asimismo, se tiene 
un efecto en las líneas de flujo, pero no es ob- 
jetivo de este trabajo analizar los límites en los 
cuales afectan estos cambios el comportamiento 
hidráulico. 


Conclusiones 


La importancia de este estudio radicó en 
conocer las pérdidas de carga presentes en 
reactores de electrocoagulación de flujo pis- 
tón a presión, de manera que se realizara su 
extrapolación para el diseño de este tipo de 
sistemas. Sin embargo, del estudio se concluye 
que se deben extender los casos de análisis, a 
fin de tener información suficiente, incluyendo 
diferentes áreas efectivas de flujo y longitudes 
de módulo. 

Con los modelos empíricos se calcularon las 
pérdidas de carga por cortante y se anexaron las 
pérdidas locales. Los errores de los resultados 
numéricos con respecto a los obtenidos en el 
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Figura 12. Promedio de errores absolutos en función de los modelos de sección anular en el R5. 
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laboratorio aumentan al tratarse de regímenes 
en transición y turbulento. 

Las pérdidas de carga consideran que el 
sistema trabajó a tubo lleno. Lo anterior es una 
simplificación del problema real debido a que 
durante la reacción electroquímica se generan 
coagulantes y gases, los cuales no se consideran 
para el cálculo y podrían impactar en la pérdida 
total. Por otro lado, las propiedades del material 
utilizado como electrodo varían con el tiempo y 
uso, esto implica que los coeficientes de fricción 
varían de igual manera, por lo que es necesario 
considerarlo en el cálculo de las pérdidas. 

Las pérdidas de carga calculadas con lo ante- 
rior se compararon con los datos recabados en el 
laboratorio, obteniendo los errores absolutos del 
ajuste de cada modelo. Como resultado inicial 
se tuvo que los modelos para sección circular 
tienen un mejor comportamiento, sin embargo, 
esto es sólo para las secciones anulares; para 
el caso de secciones anulares múltiples, son 
los modelos anulares los que tienen un mejor 
comportamiento. Luego de revisar los errores 
presentes, se concluye que la combinación de 
modelos Poiseuille-Bahrami et al. (P-BR) es la 
que ajusta a los datos de experimentación con 
un error absoluto de 0.79 m. 

Como herramienta de aproximación, el algo- 
ritmo generado se ajusta para flujos laminares y 
de transición; sin embargo, en flujos turbulen- 
tos no es aplicable el algoritmo. Por otro lado, 
dentro de las simplificaciones para el cálculo de 
los factores en los reactores de sección anular 
múltiple, se consideró que las velocidades en 
cada sección eran iguales, por ello se propone 
revisar esa simplificación por medio de modelos 
numéricos o paquetería de dinámica de fluidos 
computacional (CFD). 

Se deja para estudios posteriores ampliar el 
análisis de geometrías analizando el impacto de 
la relación sección de flujo y perímetro mojado 
en la aplicación de los factores de fricción de 
manera conjunta. 
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Variables 

Cv Coeficiente de Hazen-Williams 
Es Coeficiente de Swamee é Jain 
D Diámetro de tubería 

D, Diámetro hidráulico 

d, Diámetro interior 

d, Diámetro exterior 

ed Factor de fricción 

Por Factor de forma 

8 Aceleración gravitacional 

k Relación de diámetros 

L Longitud de tubería 

mm Coeficiente de rugosidad de Mamning 
Q Caudal 

Re Número de Reynolds 

y Viscosidad cinemática 
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